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We compare the quantum and the classical description of 
the two-dimensional motion of electrons subjected to a per- 
pendicular magnetic field and a one-dimensional lateral su- 
perlattice defined by spatially periodic magnetic and electric 
fields of large amplitudes. We explain in detail the compli- 
cated energy spectra, consisting of superimposed branches of 
strong and of weak dispersion, by the correspondence between 
the respective eigenstates and the "channeled" and "drifting" 
orbits of the classical description. 



I. INTRODUCTION 

In the last decade there has been a constant inter- 
est in the transport properties of the periodically mod- 
ulated two dimensional electron gas (2DEG). In partic- 
ular, in the presence of a lateral modulation of a one- 
dimensional character the resistivity may be strongly 
anisotropic, which essentially reflects the anisotropy of 
the electronic states. Two types of modulations can be 
achieved in the experimental devices: electrostatic poten- 
tial modulatipiteEra and, more recently, magnetic field 
modulations Jj^a Weak modulations of both types lead 
already to pronounced magnetoresistance effects in the 
presence of an average magnetic field Bq applied perpen- 
dicular to the 2DEG. These effects occur at low and in- 
termediate -Bo values, well below the magnetic quantum 
regime where Shubnikov-de Haas oscillations appear. At 
very small values of Bq a pronounced positive magne- 
toresistance is observed, followed at intermediate Bq val- 
ues by the "Weiss oscillations" due to commensurability 
effects. Both effects are adequately understood within 
a classical transport calculation based on Boltzmann's 
equation, and can be traced back to the. pxpdominancc of 
different types of classical trajectoriesJ3ilJ 

The positive magnetoresistance is understood as 
caused by "channeled orbits" which exist if the modu- 
lation is sufficiently strong or, equivalently, the average 
magnetic field is sufficiently small. For electric modula- 
tion they occur near the minima of the modulation po- 
tential ("open" orbitalSl), and for magnetic modulation 
near the lines of vanishing total magnetic field ( "snake" 
orbitsllj). They are always confined within a single pe- 
riod of the modulation, which we choose in x direction. 
They are wavy trajectories allowing for fast motion of 
electrons with velocities within small angles around the 
direction of translational invariance (y direction) . These 



channeled orbits occur in addition to the "drifting or- 
bits" , which are self-intersecting trajectories with loops 
(along each of which the direction of the velocity changes 
by 27r) , so that usually a low drift velocity in y direction 
results. For sufficiently small Bq^ drifting orbits may ex- 
tend over many periods of the modulation. At sufficiently 
large Bq (sufficiently small modulation amplitudes) only 
the drifting orbits survive. The "Weiss oscillations" man- 
ifest a commensurability effect depending on the ratio 
of the extent of drifting orbits (at the Fermi energy) 
and the modulation period. With increasing modulation 
strength, the positive magnetoresistance becomes more 
pronounced and extends to larger Bq valuea^uppressing 
progressively the Iow-Bq Weiss oscillations. E£l This effect 
is well understood within the classical calculatiorO, if 
both types of trajectories are adequately included, and it 
has recently also been obtained by a quantum calculation 
for a strong modulationt^l. 

A qualitatively new type of magnetoresistance effect 
has recently been observed by Ye et alt3 on samples 
with an extremely strong magnetic modulation. Sam- 
ples with a surface array of ferromagnetic micro-strips 
were measured in tilted magnetic fields, so that the ap- 
plied magnetic field had a large component parallel to the 
surface, producing a large magnetization of the ferromag- 
netic strips, while only the small perpendicular compo- 
nent determined the average magnetic field Bq in the 
2DEG. In this way a huge positive magnetoresistance 
with superimposed Shubnikov-de Haas like oscillations 
was obtained at low values of the average magnetic field, 
at which no magnetic ouantum effects should be expected 
for weak modulation. Efl It rather seems that the quan- 
tum oscillations are induced by the large-amplitude peri- 
odic magnetic modulation field. Such conditions require 
a quantum transport theory and, as a first step, the un- 
derstanding of the quantum electronic states of a 2DEG 
with a strong magnetic modulation. This is the motiva- 
tion of the present work. 

Channeled and drifting quantum states in linearly 
varying .rnjignetic field are already discussed by other 
authors.E3E^' 



I The Schrodinger equation for periodic mag- 
netic fields alternating in sign, has been solved previously, 
but only for the case when the average field is zeroE3. 
In the present paper we shall study the quantum elec- 
tronic states in strong periodic magnetic fields with a 
non-vanishing average, and compare it with the case of 
a strong electric modulation. In both situations, rather 
complicated energy spectra are obtained, with striking 



qualitative similarities and quantitative differences. For 
the case of strong electric modulation, such apGomplicated 
energy spectrum has recently been publishedli^, but with- 
out an attempt of an explanation. Wc will demonstrate 
in this paper that a close comparison with classical mo- 
tion leads to a detailed and intuitive understanding of 
these spectra and the corresponding eigenstates. 

In Sec. II we start with some general remarks on the re- 
lation between quantum and classical description of the 
2D electron motion in ID lateral superlattices, and we 
introduce suitable reduced units. In Sec. Ill we focus 
on the effect of a simple harmonic magnetic modulation 
of arbitrary strength. In Sec. IV we include an electric 
modulation, which requires a somewhat different analyti- 
cal procedure. The inclusion of electric modulation seems 
also necessary from the experimental point of view, since 
the ferromagnetic strips on the sample surface introduce 
a periodic stress field in the sample, which acts as an 
electric modulation on the 2DEG. Finally, in Sec. V we 
summarize the essential features derived in the paper and 
extend the discussion beyond the model of simple har- 
monic modulations. 

Some of the present resuLts. have been recently pub- 
lished in a preliminary form. Ell 



II. GENERAL REMARKS 

We consider a (non-interacting) 2DES in the x-y plane 
subjected to a magnetic field with z component Bz{x) — 
Bf) + B„i{x) and an electrostatic field in x direction lead- 
ing to a potential energy U{x). Our aim is a close com- 
parison of the classical and the quantum description of 
the electron motion (in terms of orbits and wavefunc- 
tions, respectively) in such fields, especially in the case 
that U{x) and Bm{x) are periodic in x with the same 
period a and vanishing average values. 

To evidence the translation invariance in y direction in 
the (either classical or quantum) Hamiltonian 
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we describe Bz{x) by an a;-dependent vector potential 
A(a;) = ^(a;)ey with A{x) = xBQ + Am{x) and Am{x) = 
L dx' Bm{x'). Then y is a cyclic variable and the canoni- 
cal momentum py is conserved, and one obtains an (one- 
dimensional) effective Hamiltonian H{Xq) — p'^/2m + 
V{x;Xo). For Bq ^ 0, the effective potential can be 
written 



V(x-X^)^-ujI 



Xn 



Am{x) 

Bo 



+ U{x) , (2) 



where Xq — —py/eBo is the center coordinate of the 
effective potential and uq = eJSg/m is the cyclotron fre- 
quency, both in the absence of modulation. 

In the quantum description, the reduction to a one- 
dimensional problem is achieved by the product ansatz 



1/2 

*«,Xo(a;,2/) = Ly ex.p{ipyy/h)ilJn,Xo{x) for the energy 
eigenfunctions, where Ly is a normalization length, and 
the discrete quantum number n — 0, 1, 2, . . . counts the 
nodes of the reduced wavefunction '0ri,Xo(2;)- If V{x) and 
A„i(x) are bounded, the '^n.Xoi^) drop Gaussian-like for 
|a;| -^ oo, and for a fixed value of the quasi-continuous 
quantum number Xq the energy spectrum £'„(Xo) is dis- 
crete. 

In the classical description, we use the equation mvy — 
Py~\-eA(x\ which may also be derived directly from New- 
ton's equation, to eliminate the velocity Vy. The effective 
motion in x direction is determined by H(Xq) — E. Sim- 
ilar to the wavefunctions, the orbits for given constants 
of motion, Xq and E, are bounded in x direction, how- 
ever the energy .E is a continuous variable. For a given 
E = Ep, each position x (with U{x) < Ep) is the turning 
point of twQ-orbits which are characterized by the center 
coordinatestil 
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obtained from H{Xo) = Ep for Vx = Px/fn = 0. Here 
Rq = vp/ujQ is the cyclotron radius of electrons mov- 
ing with energy Ep = mvp/2 in the magnetic field Bq. 
For given Ep and Xq, orbits exist in intervals in which 
^o{x) ^ -^0 ^ ^0^(2;) holds. This allows a convenient 
classification oiJkhe possible orbits at fixed energy Ep and 
for varying Xo.llil Of course, the same classification can 
also be done by directly investigating the effective poten- 
tial. This may be preferred if one is interested in orbits 
at different energies but the same Xq. 

The calculation of the orbits is a simple textbook prob- 
lem, but must in general be done numerically. In accor- 
dance with the translational symmetry of the problem, 
we will in the following not distinguish orbits which differ 
only by a rigid shift in y direction. If an electron is at 
time ti at position [xi,yi) on an orbit characterized by 
the constants of motion Ep and Xq, with turning points 
XI and Xr [xi < Xi < Xr), it moves towards one of the 
turning points so that at time 



t{x;Xo,Ep) = U 



dx' 



\vx{x'-Xo,Ep)\ 



(4) 



it is at position (x,y{x]XQ,Ep)), with 
y{x;XQ,Ep) = y, -l- ' "- 



Vxix';Xo,Ep) 



dx' , (5) 



whe re \vx{x;Xo,Ep)\ = vp y^ 1 ~ V {x; Xp) / Ep ^ 
uo^[X+ix) - Xo][Xo - Xoix)] and 

Vy{x;Xa) = {uJo/2)[X+{x) + X'ix) - 2Xo]. If at one 
of the turning points Xi or Xr, where Vx{x; Xq, Ep) — 0, 
the derivative dV{x; XQ)/dx vanishes, we call this turn- 
ing point and this orbit "critical". At critical turning 
points the integrals (0) and {m diverge, so that the crit- 
ical orbits there asymptotically approach straight lines 



parallel to the w axis. For non-critical orbits the inte- 
grals (H) and (^ converge as x approaches the turning 
points, and the total orbit can be composed out of right- 
running (vx > 0) and left-running (vx < 0) pieces with 
finite traverse time T{Xo,Ef) — J^'' dx/\Vx{x;Xo,Ep)\. 
The probability density of finding the electron at position 
X is Wix;Xo,Ep) = l/[T{Xo,Ep)\vxix;X„,EF)\]. This 
is the classical analog to \'4'n,Xo{x)\'^- 

If U{x) — U{x A- a) and Bm{x) — Bm{x + a) are pe- 
riodic with period a, as we will assume in the follow- 
ing, the effective potential, Eq. (Q), has the symmetry 
V{x + a; Xq + a) — V{x; Xq). As a consequence, the en- 
ergy spectrum is also periodic, En^Xg + a) = £'„(Xo), 
and can be restricted to the "first Brillouin zone" < 
Xq < a. The eigenfunctions can be taken to satisfy 
ipn,Xo+aix) — ipn,Xo{x — a). The corresponding classi- 
cal symmetry is that an orbit characterized by Ep and 
Xq + a differs from that characterized by Ep and Xq only 
by a rigid shift of amount a in a; direction. 

The dispersion of the energy bands i?„(Xo) implies a 
group velocity in y direction. 
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which is the expectation value of the velocity operator in 
the energy eigenstate ipn,Xo ■ It is the quantum equivalent 
to the classical drift velocity, i.e. the average velocity (in 
y direction) along the corresponding classical orbit. The 
drift velocity in x direction vanishes, since the orbits are 
bounded in x direction. 



A. Suitable units 

For an economic comparison of classical and quantum 
aspects it is important to use suitable length and en- 
ergy units, which are meaningful for both the quantum 
description and the classical limit. Doing so, we will see 
that the classical features depend on fewer scaled parame- 
ters than the quantum ones. To be specific but still rather 
general, we assume in the following periodic modulations 
of the form B„,{x) = B^b{Kx) and U{x) = Vqu{Kx) 
for the magnetic and the electric modulation, respec- 
tively, where b{^) and u(^) are dimensionless periodic 
functions with period 27r and vanishing average values. 
Thus Brn,{x) and U{x) have the same period a = 2t:/K, 
but may have different shapes and phases. In the numer- 
ical examples we will use for both 6(^) and u(^) simple 
cosines, eventually with a phase shift. 

The average magnetic field Bq sets, with the magnetic 
length Iq = ^/h/JjmjQ) and the cyclotron energy Twjq, 
both a length and an energy scale, which are useful for 
quantum calculations, but have no meaning for the clas- 
sical motion. For the discussion of commensurability ef- 
fects, such as the Weiss oscillations, the cyclotron orbits 
must be compared with the period a of the modulation. 



Therefore a is a natural choice for the lengths unit. The 
choice of a suitable energy unit is motivated as follows. 

Classically, Bq determines only the cyclotron frequency 
ciJo, and one needs an independent length I to define 
an energy scale Kiag = mujQp /2. Using / as length 
unit, we may define dimensionless variables ^ — x/l and 
^0 = Xq/1. The effective potential, Eq. (^), then can be 
written as V{x]Xq) = Knagi^(^;^o) with 

i)(e; Co) - K - Co + s a{KlO/Klf + w u{KlO , (7) 

where s = B^J Bq, a(C) = /^^ dC'b{C,'), and w = Vo/Knag. 
In the quantum description, the kinetic energy operator 
— {Ti^/2m)(P/dx'^ — —El (P/d^"^, introduces a new energy 
scale El — Ti^ /{2mP), which has no classical analog. In- 
troducing the energy ratio a = i?;/Knagi we write the 
effective Schrodinger equation as 



de 



^'(C;Co)-e«(Co) 



^n.Co(O=0, 



(8) 



with £„(Co)=^n(-'^o)/Knag and '0„,4o(O = ^V'n,Xo(2;)- 

If we would take I = Iq, we had Knag = Ei = hujQ/2 
and thus simply a — 1. The effective potential Eq. (M) 
would then depend on the constant of motion ^o and, 
in addition, on three dimensionless model parameters, s, 
w, and KIq. To specify an eigenstate or, in the classi- 
cal description, a trajectory, one further needs an energy 
value e as a second constant of motion. A description 
that, for fixed constants of motion, needs three parame- 
ters to specify the effective potential and, furthermore, 
relies on Iq and hcoQ, which have no meaning in classical 
mechanics, is rather clumsy and not acceptable. 

Instead we take I — 1/K and, therefore, l^nag = Kycj 
where T4yc = mujQ/{2K'^) is the energy of a classical cy- 
clotron orbit of radius l/K in the homogeneous magnetic 
field Bq. Now the effective potential Eq. M) depends 
only on the two dimensionless modulation strengths s and 
w = Vb/T4yc, which both are well defined within the clas- 
sical approach. Also the constants of motion, ^q = KXq 
and e — E/Vcyc = [KRqY , including the dimensionless 
version of Eq. (H), 



^0^ (0 = e + s a(0 ± ^e-wu{S,) , 



(9) 



remain meaningful in the classical limit. This choice of 
units will also be very useful for a systematic discussion 
of the quantum mechanical energy spectra. Quantum 
mechanics enters the effective Schrodinger equation (Ig) 
only via the parameter a = (IqK)^, which scales the ki- 
netic energy. It determines the only true quantum aspect 
of the spectrum, namely the spacin g of the energy levels 
En (Co)- We will see in Sec. [II B that all the essential 



structural features of the energy spectrum, e.g. the com- 
plicated back-folded structure due to the coexistence of 
"channeled" and "drifting" states, are determined solely 
by the "classical" parameters s and w. The density of 
the quantized levels iniCo), on the other hand, increases 
with increasing ratio u/Iq. 



As a simple example one may consider the well known 
case of a weak electric or magnetic cosine modulation, 
which -Leads— to modified Landau bands of oscillatory 
width. QQila'tj The band width assumes minima near the 
"flat band" energies E^ = m{ujoa)^{X ± f/4)/8, with 
"+" ("— ") for magnetic (electric) modulation and A = 
1, 2, .... These flat band energies are distinct multiples 
of our energy unit Kyc, and occur at if = 7r^(A ± f/4), 
independent of the special values of the model parame- 
ters Bq and a. The level spacing, on the other hand, is 
of the order hujQ and depends in our units on hoJo/Vcyc = 



III. MAGNETIC COSINE MODULATION 

We first consider a pure magnetic modulation, U{x) = 
0, Bra{x) = B'^b{Kx), so that the effective potential 
Eq. (0) becomes 



nC;eo) = K.yc[e-eo + sa(0]' 



(10) 



For s — one obtains the well known Landau levels and 
the Landau oscillator wave functions, fnXoi^)- We use 
the set fnXa as the basis of our Hilbert space in order 
to obtain numerical solutions for s 7^ 0, by numerical 
diagonalization of H{Xo). The electron effective mass is 
that of GaAs, m = 0.067too. We further assume spin 
degeneracy. For the numerical parameters chosen here 
the size of the basis will vary between f 50-300 Landau 
levels. 

Before discussing the numerical results we summarize 
some properties of the effective potential and of Eq. (0), 
which now reduces to 



^^iO=^ + -^ci{0±KRo. 



(11) 



For a fixed ^0 the local extrema of the effective potential, 
given by dV{£_, S,o)/d£_ = 0, are the points where the total 
magnetic field is zero, i. e. the roots of 



f + s6(0 = 0. 



(12) 



and the points where the effective potential is zero, i. e. 
the roots of 



e-eo + sa(0=0. 



(13) 



An important aspect for the following discussion is that 
the roots of the first kind, Eq. (p^ , if existent, are inde- 
pendent of ^0, while those of the second kind, Eq.(13), 
do depend on ^o- We will see that orbits with ^ values 
near roots of the first kind are channeled, while those 
with ^ values near roots of the second kind are drifting 
orbits. The analytic dependence of the effective poten- 
tial on the relevant position coordinate ^ is determined 
by the modulation strength s. Therefore the number of 
its possible zeroes, the classification of orbits and the en- 
ergy spectrum depend critically on the parameter s. To 
demonstrate this, we choose in the following examples 
KO = cos^, and consequently a(^) = sin^. 



A. Weak modulation, s < 1 

For s < f , the effective potential has exactly one min- 
imum of the second kind for each value ^0 , which is due 
to the confinement by the average magnetic field. The 



functions ^q (^) in equation (pTj) have no extrema. For 
each value ^0 they determine exactly one orbit, which 
is a drifting cyclotron orbit. By this we mean a self- 
intersecting orbit consisting of loops along each of which 
the azimuth angle in velocity space, ip — arctan^Vy/vx) 
increases by 27r. A typical example is illustrated in Fig. |f| 
for s = 0.5, ^0 = ''■/2 (i.e. Xq — a/4), and two energy 
values Ep. Figure ^a) shows the effective potential. For 
a given energy E = Ep (horizontal line) a classical orbit 
exists where V{^;^o) < Ep. Figure |](b) shows the loca- 
tion of the turning points as the crossing points of the 
horizontal line ^0 = tt/S with the functions <^o (<^). The 
corresponding drifting orbit exists in the interval with 
^(7(0 < So < ^iO- The orbits in real space are illus- 
trated in Fig. |f|(c). In Fig. g we plot the corresponding 
quantities for s = 0.5 and So = i" (i-e. Xq = a/2). In this 
case the effective potential is symmetric with respect to 
the center coordinate Xq. As a consequence, the orbits 
are closed and their drift velocity in y direction is zero. 

For small energies, Ep/Vcyc — (KRo)^ < tt^ (i.e. 
2Rq < a), the extents of the orbits in x direction are 
smaller than a modulation period and essentially deter- 
mined by the local values of the total magnetic field. At 
high energies, iSp/Vcyc ^ 1, the orbits extend over sev- 
eral periods of the modulation and the extent of an orbit, 
i.e., the width of the effective potential valley at the cor- 
responding energy, is determined by the cyclotron radius 
in the average magnetic field (xr — xi ~ 2Rq). 

In Fig. g(a) we display the first 50 energy bands En^^ 
calculated from the (first 150) original, degenerated Lan- 
dau levels, for s=0.5. The level spacing of the lowest 
energy bands is seen to follow the local value of the 
total magnetic field. Fig. 0(b). This is expected from 
the local approximation -En^o ~ [n + l/2)heB{^Q)/m, 
which is valid if the extent of the wavefunctions V'n.Xo i^) 
is smaller than the modulation period. With our en- 
ergy unit Kyc the apparent level spacing of energies 
which are independent of the period a becomes propor- 
tional to Y^. For example, if the local approximation 
En^g w (n + l/2)huj{^o) holds for En^„ < 41/cyc, as 
in Fig. y(a), this implies that it holds for n + 1/2 < 
4:Vcyc/[huj{^o)] = 2[ujo/uj{^o)]/yi^. Thus, the number of 
bands which are well described by the local approxima- 
tion increases quadratically with increasing modulation 
period a. 

The local approximation fails at higher energies when 
the width of the wavefunctions becomes larger than the 
period of the modulation, and the structure of the en- 
ergy spectrum changes. Indeed it is well known from 
the limit of very weak magnetic modulation, s <C f, 
that in contrast to this local approximation the bands 
become flat at the energies Ex/Vcyc — 7r^(A + 1/4), for 



A = 1, 2, . . .HtM These flat band conditions are the 
quantum equivalents to the classical commensurability 
conditions leading to the Weiss oscillations in magneto- 
transport, and do not change their positions in a plot like 
Fig. 0(a), even if we change the modulation period. A 
larger modulation period a just leads to a higher density 
of the energy bands. 

In Fig. g(c) we plot for ^o — 7r/2 the effective poten- 
tial and the square of the energy eigenfunctions for the 
energy values considered in Fig. |l|. Width and location 
of the wavefunctions in the effective potential is in close 
agreement with that of the corresponding classical or- 
bits. In Fig. H(d) we plot the corresponding quantities 
for the symmetric situation ^o = tt, to be compared with 
Fig. g. These wavefunctions belong to (relative) extrema 
of the energy bands, and thus have zero group velocity, 
in agreement with the zero drift velocity of the corre- 
sponding classical orbits. The wavefunctions in Fig. @(c) 
belong to finite energy dispersion and describe motion in 
the positive {n=3) and the negative (n=43) y direction, 
respectively, in agreement with the correponding classical 
orbits in Fig. nl For large quantum numbers n and weak 
modulation the group velocities can be shown to reduce 
quantitatively to the drift velocities of the corresponding 
classical orbits. 

In Fig. ^ we consider the "critical" situation s = 1. 
The derivatives ^Q'i^ex) = and ^Q"{^ex) — vanish for 
^ex = (2p + l)7r (p integer), i.e. for the positions where 
the magnetic field vanishes, Eq.dlJ). For all values of ^o 
the effective potential (|j) becomes flat at these points £,ex 
(see Fig. 0(c)). The classical situation is as for s < 1 with 
the exception that for ^o = £,ex i KRq there are "criti- 
cal" orbits which asymptotically approach straight lines 
parallel to the y axis on their left (for -I-) or their right 
(for — ) side, where B{x) = 0. The dashed fines plotted 
over the energy spectrum of Fig. ^(a) show the evolution 
of the flat regions of the effective potential with ^O; i-e. 
the parabolas resulting from V{£,exT^o) with p = and 
p = ±1. In the first Brillouin zone these lines are seen as 
the back-folding of the lowest parabola centered on £^ex 
with p — 0, and they are an indication of a kind of a free 
electron motion along the lines where the magnetic field 
is zero. Close to these parabolas the energy bands have 
large dispersion near inflexion points, and the energy 
separation between adjacent bands is minimum. Simi- 
lar features have been obtained in the energy .SBectra for 
single magnetic wells by Peeters and Matulis.Ey In other 
words, such states experience a weak effective magnetic 
field, due to the constant effective potential over a sub- 
stantial spatial region. The wavefunctions corresponding 
to states with large energy dispersion have large ampli- 
tudes at the positions of flat effective potential (vanishing 
total magnetic field). This is demonstrated for two se- 
lected states [(n=43, ^o=3.18) and (n=5, ^o=0.50)] in 
Fig. H(b), together with the probability distributions of 
the corresponding classical orbits. The effective poten- 
tials together with the corresponding classical orbits are 
plotted in Fig. ^(c). The trajectory corresponding to the 



state (n=5, ^o=0.50) is close to a critical orbit with a 
critical right turning point. This leads to an enhanced 
probability density near that point, which is also reflected 
in the quantum mechanical probability density. 

We will see that for slightly stronger modulation a new 
type of nearly free motion occurs with energies close to 
the parabolas V{^ex',^o) in the energy spectrum. 



B. Strong modulation, s > 1 



For s > 1, ^^'(C) = at a{0 = cos^ = -1/s, and 
^0 (0 ^^^ extrema at the following positions: 



maxima : 
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1)7T-S, 
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where p is an integer and S = arctan y/s"^ — 1. The values 
at these extrema are 



?0 \><max ) 



(2p- 
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l)7r 
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g{s) ± KRo 
g{s) ± KRo , 



where 



g{s) = v s^ — 1 — arctan v s^ — 1 > . 



(15) 



(16) 



The effective potential V{^;£,o) has extrema of the first 
kind, Eq. (|l2), at the same positions. The extrema with 

values T^(^^„; Co) = Vcyc[{2p+l)Tr- g{s) -^o]'^ are min- 
ima if {2p + l)7r > Coj Enid maxima otherwise, and those 

with values F(CmL;Co) = Kyc[(2p+l)7r + .g(s) -^o]^ are 
maxima if (2p + 1)tt > S,o, and minima otherwise. 



1. Classical approach 

The number of zeroes of ^q {£,) — ^q depends on both s 
and <^o- If g{s) < ""■, Co (C) ~ Co has at most three zeroes. 



t(p) 



t(p) 



If Co = Co (0 for any C satisfying ^!^L < C < Cm' 
i.e. ff (2p + 1)tt ~ g{s) < Co T KRq < {2p + l)7r + g{s), 
C^(C) — Co has three zeroes. The same argument holds 
for Eq. (p^, i.e. the effective potential has three zeroes. 
For (2p - 1)^ + gis) < Co T KRo < (2p + l)7r - g{s), on 
the other hand, there exists only a single zero. 

In Fig. H we show, for s = 2 [ i.e. (?(s)=0.685], an ex- 
ample where the effective potential has a single zero near 
C/27r ~ 0.1, so that for sufficiently low energy only a 
single drifting orbit exists. The number and the type of 
the possible orbits depend on the energy. At the highest 
energy shown in Fig. 5|(a) two orbits exist (solid lines in 
Fig. 0(d)). There is a drifting cyclotron orbit extending 
over more than two periods of the modulation, with the 
left turning point on Co (C) (uppermost curve in Fig. H(c)) 
near ^/2tt = —1.1, and the right turning point on Co~('?) 
(bottom curve in Fig. @(c)) near C/27r — 1.3. Near the 
relative minimum of Co^(C) close to S,/2tt — 1.7, which 



corresponds to a relative minimum of the effective po- 
tential (thick dashed line in Fig. ^(a)), there exists a 
"channeled orbit" moving in positive y direction. We 
define channeled orbits as trajectories which have both 
turning points either on ^q (^) or on £^q (^) , in contrast to 
the drifting orbits with one turning point on ■fp (^) and 
the other on £,^{£,)- I^i contrast to the self-intersecting 
drifting orbits, the channeled orbits are always confined 
to less than a single modulation period, and they move 
without self-intersections in a relatively narrow interval 
of angles around the positive or the negative y direction 
[see Fig. [^(d)] . Note that the curvature of the trajectories 
changes sign at the positions where the total magnetic 
field vanishes, see Fig. ^(b). 

If we lower the energy to E /Vcyc = 40, we arrive in 
Fig. ^ at a situation where only a single drifting orbit 
exists (dashed lines). In general, the extent in ^ direction 
of the drifting orbits decreases with decreasing energy. 
At the lowest energy indicated in Fig. || (lowest dotted 
line in (a) and innermost lines in (c)), we have again 
a drifting orbit near f/27r = 0.1 and a channeled orbit 
around ^/27r ~ 0.6. At this low energy, the extent of the 
drifting orbit is considerably smaller than the modulation 
period. 

In Fig. we show, for the same modulation strength, 
s = 2, a situation, ^o = ""j where the effective poten- 
tial has three zeroes, as is emphasized in the inset of 
Fig. ^(a). These zeroes are separated by two shallow 
maxima. If a (positive) E value below these maxima is 
chosen, one finds three narrow drifting cyclotron orbits 
located around the zeroes of the effective potential (solid 
lines). For higher energies one may find either one drift- 
ing and two channeled orbits (dotted lines) or a single 
drifting orbit (e.g. for 0.5 < E/Vcyc < 30, not indicated in 
the figure) . Actually the "drifting" orbits located around 
S. — TT have zero drift velocity due to symmetry reasons. 

In summary, for < g{s) < tt we find for given values 
of the constants of motion, ^o and E, at least one and 
at most three orbits. For larger values of s = B^^JEq, 
more orbits may exist for a given pair of ^o a-nd E values. 
A careful analysis of the extrema of the functions ^q (^) 
shows, e.g., that for tt < g{s) < 2tt between three and 
five orbits belong to the same pair of ^o and E. We will 
come back to this case below. 

Apparently the plots of the effective potential V{^; fo) 
are very useful to see which orbits are possible for a fixed 
value of the center coordinate ^o and different energies. 
Channeled orbits exist in side valleys near relative min- 
ima of V^(C;Co)- Ifi on the other hand, the energy of 
the motion is given, the plots of the locations of turning 
points Co (0 i^ very useful to classify the possible orbits 
for different values of ^o- Channeled orbits exist near 
relative minima of ^q {^) and relative maxima of ^q {£,) . 



2. Quantum calculation 

The energy spectra become more complicated in the 
case s > 1, Fig. |^ and ||. Regions of different charac- 
ter can be distinguished in these spectra. Areas, where 
the energy bands are nearly parallel lines with low dis- 
persion (region I), alternate with regions, where steep 
bands with large dispersion seem to cross bands with 
weak dispersion (region II). In fact the energy bands 
never cross each other and the apparent intersections are 
anti-crossing points with exponetially small gaps. 

The boundaries of these regions are given by classi- 
cal values only. If the energy is scaled by the classi- 
cal cyclotron energy Vcyc, for fixed modulation strength 
s the regions II are surrounded by the parabolas E — 
^(^nl;fo) and E = F(d''ix;Co) which, for | p |< 2, are 
indicated by dashed lines in the spectra. For any fixed 
Co such a pair of parabolas gives the minimum and the 
maximum value of a certain side valley of the effective po- 
tential F(C; Co) (extrema of the first kind, see Eq. ([12|)). 
The energy interval in between these values indicates the 
depth of that valley, i.e. an energy range in which clas- 
sically channeled orbits exist, in addition to the drifting 
orbits. 

In Fig. CTb) the effective potential is plotted for the 
symmetric case Co = ti" (dotted line), corresponding to 
the classical situation described in Fig. |g. Also shown 
are the states for n = (lower solid line) and for n — 20 
(upper solid line) and n = 22 (upper dashed line). Ap- 
parently, state n = 20 corresponds to a classical drift- 
ing orbit, whereas n = 22 is the symmetric superposi- 
tion of two states corresponding to channeled orbits in 
the side valleys. The latter has practically the same en- 
ergy as the corresponding antisymmetric superposition 
(n = 21), which is not shown. On the scale of Fig. 0(a), 
all states n = 20, 21, and 22 seem to have the same 
energy, E/Vcyc ~ 38. The states n = 21 and 22 are hy- 
bridizations of states belonging to the branches with high 
energy dispersion and opposite sign of the group velocity 



mujQ dCo 



(17) 



Figure R(c) shows the effective potential for the asym- 
metric case Co/27r = 0.2, corresponding to the classical 
situation described in Fig. 0. Here we show six states, the 
ground state n = located near the zero of the effective 
potential, the two "channeled" states "bound" in the po- 
tential valley around C/27r = 0.6, the extended "drifting" 
state n — 2b near E — 40T4yc, the extended state n — 44, 
and the localized "channeled" state n = 45. The energies 
of all these states are indicated in Fig.0(a). The states 
which extend over more than a period of the modulation 
belong to weakly dispersive energy bands and correspond 
to the classical drifting orbits. The states belonging to 
the energy branches with strong dispersion have large 
amplitudes in side valleys of the effective potential and 
vanish practically outside these valleys. They correspond 



to classical channeled orbits. The apparent number of 
nodes of the large-amplitude parts of these "channeled" 
states increases with energy as if they were truly bound 
states in these narrow valleys. Note, however, that the 
wave functions of channeled states still have n nodes, but 
the corresponding oscillations are not observable at the 
scale of the figure. Outside the valleys, the maxima in 
between the nodes are a few orders of magnitude smaller 
than the main peaks inside the valleys. 

The number of quantized states within a given valley of 
the effective potential depends on the average magnetic 
field and the period of the magnetic modulation, even if 
the parameters T4yc and s are fixed. In Fig.^ the period 
a, or the field Bq, is too small to have states quantized in 
the low-energy triple minimum of the effective potential 
for the symmetric case of Fig|^(b) (see also Fig.@). To 
investigate this situation, we show in FigJ^(a) a denser 
spectrum for the same modulations strength s = 2. In 
Fig.g(b) the spectrum near <^o = ""■ is enlarged. Five 
energy values are indicated, and in Fig.g(c) the corre- 
sponding (squares of the) wavefunctions are plotted for 
the states n = 0, 2, 3, 4, and 5, together with the effec- 
tive potential. The antisymmetric state n = 1, which is 
nearly degenerate with n = 2, is not shown. This demon- 
strates that all the classical features have their quantum 
analog, provided the model parameters (here a) are suit- 
ably chosen. 

For the magnetic cosine modulation, the depth of the 
valleys of the effective potential, 

l^(d^ix;eo) - ^(dl;Co)|/K.vc = 4g{s)\{2p+l)7r-^o\ , 

(18) 

increases with the energy (i.e. with \p\ for fixed ^o), and 
thus more and more channeled states appear at higher 
energies. For sufficiently high energies the strips with 
channeled states in the energy spectra may thus extend 
over the whole Brillouin zone. This will also happen for 
sufficiently large s. The energy dispersion of the chan- 
neled states depends strongly, nearly quadratically, on 
^0 according to Eq.(|l^, which expresses the nearly free 
motion of the electrons on channeled orbits in y direction. 
For s > 1 and g{s) < it, the area of the regions 
II of the spectrum increases with increasing s, and the 
area of the regions I shrinks accordingly. For g{s) = tt, 

one has V {^l^{^; £,o) = T^(Cmalc^';Co) and the correspond- 
ing parabolas coincide, leaving no room for regions I. 
If the modulation is so large that g{s) > tt, drifting 
and channeled states coexist everywhere in the spec- 
trum. In Fig. we have chosen s — 5, corresponding 
to g{s) = 3.53. Close to the edges of the Brillouin zone, 
e.g. for fo/27r — 0.016, Fig. ||(c), we can identify chan- 
neled, e.g. n = 22 and n — 16, and drifting states, e.g. 
n = 19. But now these drifting states are relatively nar- 
row, confined in local minima of the effective potential 
and not in the wide potential well centered around ^q, 
which is given by the confinement due to the average 



magnetic field. This case is already known from the dis- 
cussion of Fig. Ig. The local minimum of the effective 
potential at ^ = is a minimum of the second kind, with 
vanishing potential. Consequently, these drifting states 
are similar to weakly perturbed Landau levels with en- 
ergy gaps fLeB{^o)/m, as can be observed by a careful 
look at Fig. ^(a). In the center of the Brillouin zone say 
for ^o/27r = 0.493, Fig. ^(d), the effective potential has 
three zeroes of the second kind (see Eq.(|l^)) near ^ = ^o- 
We therefore can find similar narrow drifting states, like 
n = 5 and n = 6, but also wide drifting states at higher 
energies, like n — 17 and channeled states in local min- 
ima of first kind, like n = 7. As in the classical picture, 
the velocity of these narrow drifting states is in general 
lower than that of the channeled states. 



IV. MIXED HARMONIC MODULATIONS 

For sufficiently strong mixed electric and magnetic 
modulations one expects a similar situation as for the 
strong magnetic modulation, with a coexistence of chan- 
neled and drifting orbits, and their quantum analogs. In 
the presence of an electric modulation, we have no ex- 
plicit analytic expressions for the minima of the effective 
potential, Eq. (|7|), not even for simply harmonic modu- 
lations. Nevertheless, a qualitative understanding of the 
classical and the corresponding quantum mechanical mo- 
tion is possible. For a given constant of motion ^o the ef- 
fective potential has side valleys with possible channeled 
orbits, if dV{£,; ^o)/d£^ — has more than one solution ^. 
This is the case, if the function 



, /.^x > /,N w u'(£) 

Co(0 = e + -(0-2YT^ 



(19) 



with w — Vb/Vcyc, assumes the value S,a{£,) — S.o at more 
than one ^ value. At such ^ values the effective potential 
has extrema with the values 



V{^;^oiO)/V,yc^wuiO 



u'iO 
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l + sbiO 



(20) 



cos^. 



To be specific, we choose for the following b{£,) 
a(^) = sin^ and u{^) — cos(^ +J£s)- 

Apparently, Eqs. (|l9|) and (|20|) provide a parametric 
representation of the possible relative extrema of the ef- 
fective potential in the energy-versus-^o diagram, similar 
to the dashed lines in Figs. M and g which define the re- 
gions II where "channeled" states coexist with "drifting" 
ones. Due to the symmetries ^o(C + 27r) = £,q{(,) + l-n 
and V[i + 27r;Co(C + 27r)) = V^(^;eo(C)), it is sufficient 
to consider only one period < ^ < 27r of the parameter 
^, provided the ^o(C) values are back- folded into the first 
Brillouin zone. 

For strong magnetic modulation (s > 1) the denom- 
inators in Eqs. (19) and (|2(]| ) lead to poles. Then the 



regions of type II extend to arbitrary high energies, sim- 
ilar to the case of pure magnetic modulation. A qualita- 
tively different behavior is obtained for weak magnetic, 
but arbitrarily strong electric modulation, since then the 
denominators of Eqs. (M) and (EG) remain positive (for 
s < 1). Consequently, for given values of w, s, and ips 
the possible values oiV{^; ^o(O); Eq. (|2y), are bound and 
channeled orbits can exist only below a certain energy. 



A. Pure electric modulation 

As a particularly simple example we consider a pure 
electric cosine modulation, s — 0, ips = 0- For weak 
modulation (w < 2) the function ^o(0 of Eq. ( |l9|) has 
a unique inverse, i.e. the effective potential V{^;£,o) has 
for all values of ^o only a single extremum, namely the 
absolute minimum, and no "channeled" states should be 
expected. The energY-S|jiectra for this weak-modulation 
limit are well knownHu'Ej and will not be reproduced 
here. Apart from a phase shift, they look similar to 
Fig. |(a) but with flat bands near E/Vcyc = 7r2(A- 1/4), 
for A = 1, 2, . . .. 

The corresponding classical trajectories at sufficiently 
high energies are drifting cyclotron orbits. At very low 
energies, E < Vq = wVcyc, a peculiarity occurs, since 
then the classical trajectories are captured within a single 
valley of the electric potential, with turning points given 
by Eq. (0) in the interval fc < f < 27r — ^c (modulo 2tt) 
with ^c ~ arccos(£'/Vo) > 0. For ^q values in the interval 
^c < Co < 27r — Cc these trajectories are self-intersecting 
drifting orbits, whereas for ^q < Cc and ^o > 27r — ^c there 
exist channeled orbits with Vy > and Vy < 0, respec- 
tively. The orbit with ^o = Cc approaches the left turning 
point at C = Cc with a tangent parallel to the x axis, and 
that with Co = 27r — fc does the same at the right turning 
point C = 27r — fc- This peculiar low-energy behavior is, 
of course, not restricted to the weak modulation limit, 
but occurs always when the trajectories are captured in 
a minimum of the electric potential, i.e. for E < Vq. It 
is demonstrated in Fig. |l^(b), and will not be discussed 
further. 

If the electric modulation is strong enough, w > 2, the 
function Co(C) = C ~ (■'^/2)sinC, Eq. (119), has extrema 
at C+ — arccos(2/u)) > and C- = —£.+ (modulo 2tt) 
with values Co(C±) = -fg{w/2), where g{s) is defined by 
Eq. (|l|). Then, for |Co| < g{w/2) the equation Co(C) = Co 
has three solutions C in the interval |C| < tt, which are 
local extrema of the effective potential with values 

V(C;Co(C))/Kyc = 1 + iw/2f - [1 - {w/2) cosC]' . (21) 

In order to find in the energy spectra the regions II 
corresponding to side valleys of the effective potential, 
one may proceed as follows. One plots in the extended 
zone scheme F(C; Co(C)) versus Co(C)j starting at C = — tt, 
where Co(C) = ""■ and ^(C;Co(C)) = —Vo- With increas- 
ing C, also Co(C) and ^(C;Co(C)) increase and reach at 



C = C- their maximum values g{w/2) and Kyc(l-l-u'^/4), 
respectively. As C increases from C = C- to C = 0, Co(C) 
and l^(C;Co(C)) decrease towards the values and Vb, 
respectively. Increasing C from to tt leads to the mir- 
ror image of the described trace with respect to Co = 0: 

Co(C) - -Co(-C) and F(C;Co(C)) - ^(-C;Co(-C)). Fi- 

nally these four line segments have to be folded back into 
the "first" Brillouin zone < Co < 27r to obtain the ab- 
solute minimum of the effective potential as a function of 
Co and the boundaries of the regions II. In contrast to the 
strong magnetic modulation, these regions become nar- 
rower with increasing energy and end at Co = '^9{w/2) 
(modulo 27r) with energy E/Vcyc = 1-1- w^/4. For 
g{w/2) > IT the back-folding leads to an overlap of differ- 
ent branches of the region II, that is the coexistence of 
back and forth running "channeled" states with "drift- 
ing" states in the same area of the i?-Co diagram. Fig- 
ure 10 shows for a typical example the quantum mechan- 
ical energy spectrum together with the boundaries of re- 
gion II obtained in this manner. The "very complicated! 
energy spectrum obtained recently by Shi and Szetotj 
for strong electric modulation is thus explained by the 
coexistence of chaipeled and drifting states. 

In previous worklijiid it was pointed out that, for given 
modulation period a and strength Vq and given energy 
E — Ep, channeled orbits can exist only if the magnetic 
field Bq is smaller than a critical field Bcrit- Solving 
Ep/Vcyc = I + w'^/4: for Ep > Va = wVcyc and w > 2 
with respect to the magnetic field, one obtains the known 
resultM 
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(22) 



B. Weak magnetic modulation 

If a magnetic modulation is added to an electric one, 
very complicated interference effects may result. Only 
if the phase shift ips is zero or tt, the resulting energy 
spectrum will be symmetric in Co- Even in that case, the 
distribution of channeled states (regions II) in the -E-Co 
diagram may become rather complicated, especially at 
low energies. For the mixed case channeled states may 
occur even if the modulation parameters w and s are not 
large enough to produce them for the pure electric and 
the pure magnetic modulation of these strengths. For 
weak magnetic modulation, < s < 1, and arbitrary 
strength of the electric modulation, channeled states can 
exist only below a certain energy, as in the pure electric 
modulation case. 

For arbitrary phase shift ips the energy spectrum may 
be so asymmetric that in a certain energy range only 
channeled orbits exist which carry current in one (say the 
positive y-) direction, but no channeled orbits carrying 
current in the opposite direction. Such a situation is 



presented in Fig. |ll|. The regions II, where channeled and 
drifting states coexist, is again calculated from Eqs. (|l9| ) 
and (pO|), i.e. from purely classical arguments. 

If a 2DEG is subjected to such an asymmetric mixed 
modulation, it may happen that in the thermal equi- 
librium the channeled states carry a finite current. Of 
course, this current must be compensated by a corre- 
sponding opposite current carried by the drifting states. 

For large magnetic modulation, s > 1, and arbitrary 
electric modulation, the magnetic modulation dominates 
the energy spectra at large energies. The regions II with 
channeled states become more and more important, as 
can be seen from the pole structure of Eqs. (Il^) and (20). 
For very weak electric modulation, the results reduce to 
those of the pure magnetic modulation, apart from some 
peculiarities at very low energies, where additional re- 
gions of channeled orbits may exist. 

In magnetically modulated systems prepared by de- 
position of magnetic micro-strips there is always an in- 
duced electric modulation due to the interface stress be- 
tween the ferromagnets and the semiconductor.Q The 
phase shift with respect to the magnetic modula.tion oc- 
curs when the external magnetic field is tilted El Never- 
theless, in the known experimental situations the stress 
potential amplitude is presumably much weaker than our 
bare potential. 



V. SUMMARY AND DISCUSSION 

We have discussed in detail the quantum electronic 
states and energy spectra En{Xo) of a 2DEG in strong 
one-dimensional magnetic and electric superlattices, and 
in a non- vanishing average external magnetic field. By 
comparing the quantum results with the corresponding 
characteristics of the classical motion, we achieved a de- 
tailed and intuitive understanding of the energy spectra 
and eigenstates. We found that the complicated parts of 
the energy spectra ("regions IF), where branches with 
strong dispersion coexist with those of low dispersion, 
coincide with the areas in the E — Xq diagram in which 
classically "channeled" orbits exist. 

For a systematic investigation of the possible energy 
spectra and eigen states, and of the corresponding types 
of classical trajectories, it is useful to exploit the scaling 
properties of the Hamiltonian. Then it is not necessary 
to vary independently all the basic model paramaters, 
i.e. the strengths B^ and Vq of magnetic and electric 
modulation, the modulation period a, and the average 
magnetic field Bq. If one uses suitable units for energy 
and length, V^yc and a/27r, respectively, one obtains the 
same classical results and the same gross features of the 
energy spectra (the same position of the regions II), if 
one changes the four parameters i?^, Vb, a, and Bq in 
such a manner that the two reduced modulation strengths 
s — B'^/ Bq and w — Vo/T4yc remain constant. For dif- 
ferent parameter sets with the same values of s and w, 



only the density of the energy bands is different in the 
plot oi En{Xo)/Vcyc versus KXq, not its overall appear- 
ance. This is illustrated by Figs. R(a) and @(a), for which 
the regions II coincide. The reason for this behaviour 
is that, in these energy and length units, the effective 
potential is invariant under the scaling transformation 
B?n -^ iB^m, Bo -> 7^0, a^ Xa and Vo -^ t'A^Fq, for 
arbitrary positive 7 and A. To leave the quantum result 
exactly unchanged under a change of the four model pa- 
rameters, one has to keep a = (loK)"^ also unchanged. 
This is because only with the restriction A = 1/y^ the 
kinetic energy operator is also independent of the scaling 
parameter 7 [see Eq. (0)]. Thus, in the suitable units, 
the exact quantum result depends only on three indepen- 
dent parameters instead of four, and the characteristic 
classical features depend only on two. 

There is a close correspondence between the quantum 
states belonging to strong-dispersion branches of the en- 
ergy spectrum and the classical channeled orbits. These 
orbits occur near lines of vanishing total magnetic field or 
near minima of the electric modulation potential and are 
restricted to individual side valleys of the effective poten- 
tial. They are always restricted to a part of a single mod- 
ulation period in x direction and represent a fast motion 
along (wavy) lines without self-intersections in positive or 
negative y direction. The corresponding quantum states 
are also essentially confined to the same space region and 
belong to energy branches with strong dispersion. At a 
given value of the constant of motion Xq, channeled or- 
bits exist in energy intervals bounded by adjacent relative 
minima and maxima of the effective potential, defining 
bottom and top of the corresponding side valley. Plotting 
these classically defined extrema versus Xq, one obtains 
the boundaries of the regions II of the quantum energy 
spectrum. Classically, for each channeled orbit there ex- 
ists a "drifting" orbit with the same constants of mo- 
tion Xq and E. These drifting orbits are self-intersecting 
trajectories which, for sufficiently large energy, extend 
over more than one modulation period in x direction and 
drift slowly in y direction. The corresponding quantum 
states belong to low-dispersion branches of the energy 
spectrum. Quantum mechanically, the channeled states 
do not appear at exactly the same energies as the drift- 
ing states, and they usually have a larger energy spacing 
than the latter, since they are confined to a narrower 
effective potential well. 

We have demonstrated these features by model cal- 
culations based on simple harmonic modulation fields. 
Qualitatively the obtained results and the methods to 
derive them can easily be extended to more general mod- 
ulation fields, containing higher harmonics. This will be 
necessary, if the distance of the 2DEG from the sam- 
ple surface is not much larger than the_4period of the 
surface structure creating the modulation.til Anharmonic 
effective modulation potentials may also result from non- 
linear screening effects,, even if the bare modulation po- 
tential is harmonic.E3E3 

We have also performed several additional calculations 



and consistency checks which are not documented in the 
main text. E.g., we have checked the equivalence of clas- 
sical drift velocity and quantum group velocity beyond 
the analytically accessible case of very weak modulation 
fields. For some examples with strong modulation, we 
evaluated the quantum mechanical group velocity along 
several energy bands En{Xo) and compared the result 
with the drift velocity of the corresponding classical tra- 
jectories with the same energy and Xq values. For most 
parts of the bands the two velocities agreed perfectly. 
A systematic deviation was observed only in parameter 
regimes where the classical trajectories are close to crit- 
ical orbits, which have no quantum analog. Near the 
critical orbits the modulus of the classical drift velocity 
increases rather rapidly, whereas the quantum mechani- 
cal group velocity shows no anomaly. 

We have also extended the band structure calcula- 
tions to very strong magnetic modulation {B^^/Bq = 20). 
While at high energies a complicated superpositon of 
bands with steep and with flat dispersions, similar to that 
in Fig. ^a), was obtained, the bands at low energies tend 
to cluster into groups separated by relatively large gaps. 
The low-energy part of the spectrum was already remi- 
niscent of the spectrum for vanishing average magnetic 
field, where one obtains a one-dimensional—Bloch energy 
spectrum for each value of py = —eBoXollS 

Concerning previous and forthcoming transport calcu- 
lations, we conclude from the close correspondence of the 
quantum and the classical approach, that at weak average 
magnetic fields the classical calculations are appropriate, 
provided the modulation fields are not too strong. For 
the very strong magnetic modulation mentioned in the 
introduction, it may however happen, that the energy 
level spacing of "channeled orbits" exceeds the thermal 
energy ksT in a regime where hujo <C fc^T. Then we 
would expect modulation induced quantum effects in the 
positive-magnetoresistance regime at low Bq. 
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FIG. 1. (a): Effective potential V^(5;^o) for magnetic co- 
sine modulation with s — B^/ Bq — 0.5 and S,o/2-k = 1/4. For 
a given energy E-e/Vcyc ~ (KRo)^ (horizontal lines) classical 
orbits exist where V{^;^o) < ^f- Solid line Ef = 17.6T4yc, 
dotted line Ep — 1.1814yc. (b): Locations of turning 
points 5(, (^) for the Ep values indicated in (a), same cod- 
ing. Orbits with energy Ef and ^o exist in an interval with 
^0 iO !E: Co < Co^(C)- (c): Corresponding orbits in xy-space, 
three cycles are shown each, the sense of motion is from filled 
to open dot. 



FIG. 2. As Fig. |i|, but ^o/27r = 1/2. The effective po- 
tential is symmetric and therefore the guiding center of these 
drifting orbits does not drift in (c). 



^ D. Weiss, K. v. Klitzing, K. Ploog, and G. Weimann, Eu- 
rophys. Lett. 8, 179 (1989), see also in High Magnetic 
Fields in Semiconductor Physics II, edited by G. Landwehr, 
Springer Series in Solid-State Sciences Vol. 87 (Springer- 
Verlag, Berlin 1989), p. 357. 
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FIG. 3. (a) Landau bands for s = 0.5. Bq = 0.2 T and 
a = 800 nm, so that Kyc = 0.85 meV and a = 0.041 and 
(b) total magnetic field B{^q). The marked points on Landau 
bands 43 and 3 are the states for which the wave functions 
are shown in (c) and (d) in arbitrary units together with the 
corresponding effective potentials (dashed hue). The wave 
functions are plotted with an offset, indicating the energy of 
the state. The states of (c) and (d) are to be compared with 
the classical orbits in Figs. Q and nl respectively. 



FIG. 9. (a) Energy spectrum for s = 5, -Bo = 0.1 T and 
a = 800 nm, so that Fcyc = 0.213 meV and a = 0.16. (b) 
Total magnetic field, (c) Effective potential and specific states 
for ^o/27r = 0.016. Typical channeled states n = 16 and 
n = 22 and a narrow drifting states, n = 19. (d) Effective 
potential and specific states for ^o/2tt — 0.493. Wide drifting 
state, n — 17, narrow drifting states, n = 5 and n = 6, and a 
channeled state, n — 7. 



FIG. 4. (a) Energy spectrum for s 



Bo 



0.1 T 



and a = 800 nm, so that Vcyc ~ 0.21 meV and a — 0.17. 
The dashed lines show V(Co, (2p + l)7r) with | p |< 1. (b) 
Quantum mechanical (thick lines) and corresponding classical 
(thin lines) probability densities for two states. The chosen 
states are marked with dots in (a), (c) Efi'ective potentials 



and classical orbits for Ef — TV^yc , ^o - 
Ef = 39.61/cyc,.^o = 3.18 (dotted lines), 
indicate the Fermi energy. 



: 0.5 (solid lines) and 
The horizontal lines 



FIG. 5. (a): Efi'ective potential V(^;^o) for magnetic co- 
sine modulation with s — B'^/Bq = 2, and ^o/27r = 0.2. For 
given energy Ef (horizontal lines) classical orbits exist where 
^(^iCo) < Ep. (b): Total magnetic field, (c): Locations of 
turning points ^g (5) for the Ef values indicated in (a). The 
outermost pair of lines belongs to the largest, the innermost 
pair to the smallest Ef value. The constant of motion ^o 
appears as a horizontal line in this plot (^o/27r = 0.2 is indi- 
cated). Orbits with fixed energy (i. e. fixed curves Co (0) ^^id 
this value of Co exist in intervals with Co'(C) ^ Co < Co'(C)- Or- 
bits, plotted in (d), with one turning point on ^q {^) and the 
other on Co'(C) ^re drifting orbits, the others are channeled 
orbits (see text). 



FIG. 10. Pure electric modulation of strength w = 20. (a) 
Energy spectrum (first 80 bands) for Bq — 0.05 T and a = 800 
nm, so that 'Kyc = 0.053 meV and a = 0.67. (b) Drifting orbit 
(solid line) and channeled orbit (dotted line) for Co/27r = 0.3. 
The energies (see marked states in (a)) are chosen in such a 
way, that Cc < Co for the drifting and Cc > Co for the channeled 
orbit, see text. 



FIG. 11. (a) The first 175 Landau bands for combined 
magnetic and electric modulations. Bo = 0.1 T and a = 800 



nm, so that Vex 



0.213 and a — 0.16. The modulation 



strengths are s — 0.8 and w = 14.3, the relative phaseshift 



is 7r/2. (b) Co "Curves for -B/Kj 



110. At this energy Co 



has local extrema, while there are none in Co 
channeled orbits have negative velocities Vy. 



Therefore all 



FIG. 6. As in Fig. | but for s = 2, and Co/27r = 0.5. Hor- 
izontal lines in (a) are for -Bp/Kyc = 40 and Sp/Kyc = 0.3. 
The inset shows V{^; Co) enlarged between C = and C = 27r, 
where it has three zeroes, (c): For both indicated energies 
there exist three orbits, one drifting and two channeled orbits 
for Ef/Vcjc = 40, and three drifting orbits for Ef/Vcjc ~ 0.3, 
plotted in (d). Since the effective potential is symmetric, there 
is no guiding center drift for the central drifting orbits. 



FIG. 7. (a) Energy spectrum (first 75 bands) for s = 2. 
Bo = 0.05 T and a = 800 nm, so that Kyc = 0.053 meV 
and a = 0.67. Effective potential and specific states (b) 
(n=0,20,22) for Co/27r = 0.5 and (c) (n=0,l,2,25,44,45) for 
Co/27r = 0.2. 



Vcyc = 0.851 



FIG. 8. (a) Energy spectrum (first 100 bands) for s 
Bo = 0.2 T and a = 800 nm, so that 
meV and a — 0.041. (c) Effective potential and specific 
states (n=0,2,3,4,5), marked in extract of spectrum (b), for 
Co/27r = 0.5. 
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